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We search for self-gravitating oscillating field lumps (pulsons) in the scalar model with logarithmic 
potential. With the use of a Krylov-Bogoliubov-type asymptotic expansion in the gravitational 
constant, the pulson solutions of the Einstein-Klein- Gordon system are obtained in the Schwarzschild 
coordinates. They are expressed in terms of solutions of the singular Hill's equation. The masses 
of the obtained pulsons are calculated. The initial conditions are found under which the pulson 
solutions become periodic. These conditions are then used in direct numerical integration of the 
Einstein-Klein-Gordon system. It is shown that they do evolve into a very long-lived periodic 
pulson. Stability of the self-gravitating pulsons and their possible astrophysical applications are 
briefly discussed. 
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I. INTRODUCTION 

A large number of modern astrophysical observations 
suggest the existence of scalar fields in our Universe as 
possible candidates for dark matter. Pulsons are local- 
ized configurations of the fields having oscillating energy 
density. Numerical simulations of Scidel and Suen [l| 
have revealed the existence of long-lived self-gravitating 
pulsons, so-called oscillating soliton stars or oscillatons, 
in the Einstein-Klein-Gordon (EKG) system 
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with the potential U(4>) — (m 2 /2)(f> 2 corresponding to a 
free massive scalar field. The authors have established 
that soliton stars can be formed from rather general ini- 
tial field distributions due to specific relaxation process, 
the gravitational cooling. 

Pulsons were first observed numerically by Bogolubsky 
and Makhankov 0| in the Klein-Gordon (KG) model with 
4> 4 and sG potentials. In these cases, in the absence of 
gravity, the formation of the pulsons occurs solely due to 
self-coupling effects. In the present-day literature such 
configurations are often called oscillons, but below we 
shall use their original name, pulsons [3[. 

Subsequent investigations have shown that pulsons ex- 
ist in various models and spatial dimensions, and that 
they evolve from the diversity of initial conditions [3-[2^| 
(see [13] for a review) . It turns out that pulsons can arise 
from both uniform and non-uniform field distributions. 
Thus pulsons can emerge in scalar condensates due to 
the parametric instability of the spatially uniform back- 
ground oscillating near a vacuum value EM EM HH . 
In this case the energy of the background oscillations 



is transferred to an incipient pulson via the resonance 
mechanism. Quite a different scenario is realized when 
pulsons are formed from localized field distributions that 
appear, e.g., in shrinking cylindrical domain walls [7[, in 
collapsing spherical bubbles [MEH, or at bubble collisions 
!l5i]. In such a case an initial field lump sheds excessive 
energy by radiation of scalar waves (gravitational cool- 
ing of the soliton stars) and settles into a quasi-stable 
state, the pulson, whose lifetime depends strongly on the 
initial conditions. This suggests the existence of such 
initial conditions that evolve into very long-lived quasi- 
periodic, or even infinitely long-lived periodic pulsons. 
The latter would imply the existence of exact localized 
time-periodic solutions. For the 4 , (f> 3 — 4 , and sG 
models, certain of these initial co nfig urations have been 
found numerically in [To|, El> EH, EH- Recently, in Ref. 
}24j small amplitude pulson solutions of the EKG sys- 
tem have been obtained for the potentials expansible in a 
power series. This brings up the following question: How 
does gravity affect the dynamics of the finite amplitude 
pulsons? For example, could gravity turn non-periodic 
pulson solutions into periodic ones? Consideration of fi- 
nite amplitude pulsons takes on great significance in the 
case where a scalar field potential is not expansible in a 
power series in the small amplitude limit. 

In this paper we search for pulsons in the EKG system 
(fTJ) with the potential 
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where <p is a real scalar field, m is a bare mass (in units 
h = c = 1), and a is a characteristic amplitude of the 
field which is assumed to be finite, but not too large, so 
that yea 2 <C 1, where x is the gravitational constant. 

The nonlinear KG equation with the logarithmic po- 
tential @ was first considered in quantum field theory 
by Rosen (2f| and later by Bialynicki-Birula and Myciel- 
ski [26[ . In general, for the nonlinear KG equation the 
potential @ is the only one which permits real solutions 
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of the form (j) — a(t)u(r) to exist yj. Such singular po- 
tentials currently appear in inflationary cosmology [27] 
and in some supersymmetric extensions of the standard 
model (flat direction potentials in the gravity mediated 
supersymmetric breaking scenario) [28j j . The logarithmic 
term in parentheses arises due to quantum corrections to 
the bare inflaton mass. 

The paper is organized as follows. In Sec. UH using 
the smallness of the gravitational constant, we obtain 
the approximate solution of the EKG system (p} which 
describes time-periodic pulsons of a finite amplitude in 
the Schwarzschild metric ds 2 — Bdt 2 — Adr 2 — r 2 {d-d 2 + 
sin 2 i? dip 2 ). In Sec. Mil we use the obtained solution to 
find the initial conditions for direct numerical integration 
of the system. We show that these initial conditions do 
evolve into a very long-lived periodic pulson. Stability of 
the self-gravitating pulsons and their possible astrophys- 
ical meaning are briefly discussed in Sec. IIV1 



II. SOLUTION 

After the scaling t —> t/m, r —> r/m, (fi/a — > 0, 
xa 2 /2 — > x, the system (fT]) takes the form 

^ + ^ = -r^ 2 + ^ + J 40 2 (l-ln0 2 )V (3) 
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where x <C 1 is the rescaled gravitational constant. 
Looking for localized solutions, we impose the bound- 
ary conditions <p(t, oo) = 0, A(t, oo) = 1, B(t, oo) = 
1, r (t s O) =0, A(t,0) = 1. 

If we set x = 0, from ©-(O we immediately obtain 
A = B = 1 and arrive at the nonlinear Klein-Gordon 
equation 



4>tt - 4>rr - (2/r)4> r - In 2 = 0. 
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This equation has a whole family of exact pulson solu- 
tions [3, H @] ■ The simplest of them is given by 

r)=a(*)e( 3 -'' 2 )/ 2 , (7) 

where a(t) satisfies the equation of a nonlinear oscillator, 

a tt = -dV{a)/da, V{a) = (a 2 /2)(l - lna 2 ). (8) 

As is clear from the shape of the potential V(a) depicted 
in Fig. [TJ oscillations are possible in the range —1 < 
a(t) < 1, so we shall consider below that the pulson's 
amplitude may be finite, \<p\ < 0(1). 




FIG. 1. The shape of the potential V(a). Initial conditions 
for the nonlinear oscillator (|8]) are a(0) = a ma x, at(0) = 0. 



For small x <C 1 we construct the Krylov-Bogoliubov- 
type asymptotic expansion (see, e.g., [291 ] ) near the non- 
gravitating pulson, 

0(t, r) = [a{0) + xQ(8, r) + 0(x 2 )] e^^ 2 , (9) 
e t = l + xQ + 0(x 2 ), (10) 

where a(8) satisfies Eqs. ©, with the phase 9 instead of 
t, and the initial conditions a(0) = a max < 1, ag{0) = 0. 
The function Q(8,r) and the constant fi to be found 
describe the deviation of the pulson's shape from the 
Gaussian one and the frequency shift Soj/io — xfl due 
to gravitational effects. 
Setting in Eqs. 0} 

A(t,r) = (l-r g /r)- 1 . / B(t,r) = (I — r g /r)e~ s (II) 
and using (HJ), (fT0|) . we find 

r g (t,r) = x 

Jo \ 
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s(t,r) = 2x 
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+0(x 2 ), (13) 

where a = a(6(t)), V ma , x = V(a ma _ x ). Substituting © 
into ([5]) leads to the equation for Q(6,r): 

Q88-Qrr+C2/r)(r 2 -l)Q r -(2+\na 2 )Q = S(a,r), (14) 

where 

S(a, r) = a{V max [V^(2 - r 2 - In a 2 )(2r)- 1 e r2 erf r 
+ 3r 2 - 4 - 3 In a 2 ] - 5aV + 2aV 

+ a 2 - 2a 2 In 2 a 2 }e 3 - r2 - 217a In a 2 . (15) 
Its solution is given by 



Q{6,r) = -J2c n X n (e)H 2n+1 (r), 



(16) 
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where c„ = tt" 1 / 4 [2 2n+1 (2n + 1)!] 1/2 , and H 2n+1 (r) 
are Hermite polynomials. The functions X n {9) must sat- 
isfy the non-homogeneous singular Hill's equation 

X n ge + (E- 2 -In a 2 ) X n = f n (a), (17) 

where E = E n = An, 



f n (a)=2c n S(a,r)H 2n+l (r)e' r rdr. (18) 
Jo 



The calculation gives 

/o(a) = D (a) - \/27r 1/4 rJa In a 2 , 
f n {a) = D n (a) (n = l,2,...), 
(-l)"(2n)!c„ /?F 3 
^( a )= 2n+4(2n _ 1)n! V2 e 

x{a 3 (4n 2 - l)(4n 2 -4n- 7 
-2F max a[247i 3 + 20n 2 - 46n 
+4(2n- l)(6n + 5)hia 2 ]}. 



(19) 
(20) 
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Note that f n (a) is a T-periodic function of 9, while In a 2 
on the left hand side of Eq. (|17|) is a T/2-periodic one, 
where T is a period of a{9). 

Solutions of the homogeneous singular Hill's equation 
were investigated in Ref. 1211 . In accordance with the 
Floquet theory (see, e.g., [30] ) Eq. (|T7| with f n = has 
two linearly independent solutions of the form ip(9)e^ 
and ip(—9)e~ fJ ' s , where [i is a characteristic exponent, and 
ip(6) is a T-periodic (T/2-periodic or T/2-antipcriodic) 
function. Obviously, we can set <p(0) = 1. Let X ± (9) be 
two solutions of the homogenious Eq. (fT7| (with f n — 
0) satisfying the conditions X+ (0) = 1, X+ (0) = 0, 
X" (0) = 0, X~ (0) = 1. They can be written as 



X+ (9) = 
X- (9) = 



\p{9)e^ 
1 



-pen 



p(-0)e" 



(22) 
.(23) 



If \X + (T/2)\ > 1, we have the resonance case: /x > and 
is determined by the equation cosh(/^T/2) = \X + (T/2)|, 
(p(9) is a real T/2-periodic or T/2-antiperiodic function, 
and hence oscillations of X^ (9) grow exponentially with 
9. If \X + (T/2)\ < 1, we have the non-resonance case: 
fi = iv, and ip(0) is a complex T/2-periodic function 
such that <p*{9) = <p(—0). Hence the solutions X ± (9) 
are bounded. They can be periodic (with some period), 
or non-periodic depending on v, which is determined by 
cos(i/T/2) = X + (T/2). These cases are realized in dif- 
ferent domains of the (E, a^ax) plane that make up a 
stability-instability chart. The domains with /i > are 
known as resonance zones. The special case \X + (T/2)\ = 
1 is realized on their boundaries where /i = 0. Then one 
of the solutions, either X + (9) or X~ (9), is a T-periodic 
(T/2-periodic or T/2-antiperiodic) function, and another 
one is proportional to the product of this function times 
9 plus some T-periodic function (T/2-periodic or T/2- 
antiperiodic, respectively). 



The surface fJ,(E, a^^) over the resonance zones has 
been constructed in Ref. [2lJ. For discrete E = An the 
above functions acquire the subscript n, so we shall write 

<f n (0), X± (9), //„(<«)■ 

Each cross-section of the surface fi(E, a 2 ^) with the 
plane E — An, n = 1,2,..., gives the characteristic ex- 
ponent n n as a function of a^ax- This function is repre- 
sented by a series of peaks fj, n > separated by intervals 
of stability. By superposing the curves Mn(a 2 nax ) f° r a ^ 
considered modes n = l,2,...,N, one gets the pattern 
shown in Fig. [2j The mode n = corresponds to the 
above special case \i = and thus does not contribute to 
the pattern. The obtained composite plot gives an idea 
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FIG. 2. A collection of the resonance peaks obtained by su- 
perposition of the functions fJ- n {a.max)- 

of the existence of unstable and (quasi) stable modes in 
different regions of the a^ax axis and demonstrates the 
tendency to progressively fill the interval < a 2 ^ < 1 
by the resonant peaks as the successively higher energy 
levels E n = An are accounted for. 

In terms of X^ (9) the general solution of Eq. (fT71) is 
written as 

X n (9) = (x n (0) - J x-f n de) X+ (9) 



+ \X ng (0) + X+f n d9)X-(9). (24) 



The solutions X^ 1 (9) have the form 

K {o) = e<r 1/3 (e 2 ) - & f m 2 ) de, (25) 

Jo 

A7 (9) = -{ujI - 1)-%, (26) 

where the notations £(6) = a/a max , a; 2 = 1 — lna 2 nax are 
introduced, 



e 2 = K 2 -i)(i-a<r 2/3 (£ 2 ) 



(27) 



is the first integral of Eq. ([5} in terms of £(#), and the 
functions <?(£ 2 ) and K(^ 2 ) are 



' » v wg + (l - € 2 )~ 1 € 2 in 
K(?) = ^yrfr4- 1/8 (f 2 )- 



3/2 



(28) 
(29) 
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Note that < (1 - uo 2 f' 2 < o(£ 2 ) < 1, dq/d£ 2 > 
(£ 2 < 1), and q(l) = 1. Since AT(£ 2 ) is a sign-definite 
periodic function of 9, its average K = T~ x J* Xd(9 ^ 0, 
so the solution (|25l) can be represented in the form 
x o ( e ) = -K£ 9 9+iP(6) : where V>(6>) is a T/2-antiperiodic 
function with ip(9) =0 [here and elsewhere the bar means 
the average over the period T of a(9)]. Thus oscillations 
of Xq (9) grow linearly with 9 for any a max . This is in 
agreement with the fact that Xq (T/2) = —1 and the 
line E = is the boundary of a resonance zone on the 
(£,< ax ) plane [p. The equality Jf^ (T/2) = -1 im- 
mediately follows from Eq. (|25l) if one takes into account 
that £(T/2) = -1, fr(T/2) = (see Fig. [J). 

The requirement of boundedness of the general solu- 
tion Xo (9) (|24l) determines f2 and, hence, the frequency 
shift in accordance with Eq. (fTU|) . Indeed, substituting 
X^ (0), X^ (9), and / (a) into Eq. ([24j) and integrating 
by parts, we find that the linearly growing terms cancel 
out if 



n 



V2 X+(9)D (a) 



ln< 



X (0) X 



(30) 



Under the condition (|30|) the solution Xo (0) is a bounded 
T-periodic function. 

To obtain the corresponding conditions for n ^ 1, we 
substitute (f22|) . (|23j) into (J24J) and require that X„ (6>) = 
X„(6> + T). In this equality the integrals between the 
limits and 9 cancel out. The remaining terms make 
up a linear combination of the independent solutions 
(p n (9)e^ ne and <p n (— 6>)e~ AI ™ e . Equating to zero coeffi- 
cients of these solutions and using the identities 

e-^ T ' 2 Vn {9)e^f n {a) 



^ T/2 M-8)e-^ 8 fn(a) 

(31) 



X+ (9) f n {a) 
cosh 0„T/2) ' 



we arrive at the conditions 



X„(0) 
X nfj (0)=0. 



T 



Xitn (T) 



X+ 8 (T) 
sinh (/J, n T) ' 



X+(0)/ n (a), 



(32) 

(33) 
(34) 



Note that X+ (T) 7^ because \i n + <f ng {0) 
in (|32p is proportional to the Wronskian 
(<^„((9)e^ e , ^„(-6»)e-^" e ). Equation (02) can 
be easily derived if one expresses <p n (9) from (|22[) in 
terms of X+ (0) and X+ (0 + T) and takes into account 
that X+ (T) = cosh(^„T). 

Interestingly, Eq. ([33)) is still valid on the bound- 
aries of resonance zones, Eq. (|34[) being no longer nec- 
essary. In particular, this is true for n = 0. Indeed, dif- 
ferentiation of (US]) gives X+ (T) = -TXlna 2 nax . To 



a max lna 2 lax X 0e (9) in ([19]) and, integrating by parts, 
take into account that W (x£ (9) , X^ (6>)) =1. As a 



result, we arrive at the condition (|30[) again. 

Thus, under the conditions (|3"0"[) . (j3"3"[) . (fMl) the solution 
(IT51) is T-periodic with respect to 9. This means the 
solution (O is also periodic [with the period (1 + >&i)~ 1 T 
with respect to t\. Note it involves the free parameters 
a m ax, X (0), and X Og (0). 

To be certain that the obtained solution is correct, we 
examine the mass conservation law. The mass of a self- 
gravitating field lump is defined as M — Att J q T^r 2 dr, 
where T® is the energy density of the scalar field involved 
in the EKG system ([IJ . In terms of the rescaled variables 
it can be written as M = (2k<j 2 /m) lim r ._ ) . 00 (r ff (£, r)/x), 
where r g (t,r) is defined in (IT2l . x being the rescaled 
gravitational constant. This limit must be time indepen- 
dent. To check this, we substitute the solution (JSJ) into 
(fT2| and calculate the limit oir g /x in the first order in x 
using the orthogonality of the Hcrmite polynomials. The 
result is given by 



M = 



\f2a-a 



'77-1/41/ 



X o (0)lnai 



2 

max 



128 



;i + 141na 2 nax ) 



0(x 2 



(35) 



calculate X$ (9) fo(a) we substitute aim 



which is evidently constant. 

Since x <C 1, the gravitational field created by this 
mass is weak, as is clearly seen from (fTTj) - (Tl"3")) . In the 
limit n — > the gravity vanishes. However, the rescaled 
scalar field persists, satisfying Eq. (J5J), and its ampli- 
tude may have any value in the range < a max < 1. As 
flmax changes from unity to zero, the pulson's frequency 
changes from zero to infinity, correspondingly. In partic- 
ular, in the small amplitude limit the pulson's frequency 
is w = (1 - In a^) 1/2 . 

Thus, we have obtained a three-parametric family of 
the spatially localized time-periodic solutions (|9"1)- (fT3"|) of 
the system d3])- d5j) , wherein only the smallness of the 
rescaled gravitational constant x has been. Note that 
the smallness of the pulson's amplitude, a max <c 1, is not 
assumed in the above consideration. To our knowledge, 
this is the first example of the pulson solutions of the 
EKG system, that have an arbitrary frequency. We have 
named them gravipulsons. 



III. NUMERICAL SIMULATION 

Our solution, however, is an approximate one. It was 
obtained in the first order in the gravitational constant. 
Hence its deviation from an exact solution increases in 
time, as happens with any asymptotic solution in the 
theory of nonlinear oscillations [2{|. But we can go back 
in time and take the initial state of the obtained solution 
as initial conditions for direct numerical integration of 
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the starting EKG system. As a result, we have a three- 
parametric family of the initial conditions: 

0(0, r) = a ma ^ 3 - r2)/2 + xG{r: a max ,X (0)), (36) 
&(0,r) = 2^c Xo e (0)e( 3 " r2 )/ 2 . (37) 

The function G(r; a max , X (0)) = Q(0, r) e ( 3 - r2 )/ 2 de- 
scribes admissible deformations of the initial pulson's 
profile which evolve into periodic solutions. In calculat- 
ing G we assume that a^ax belongs to one of the intervals 
of quasi-stability [HHl where X± (6>), with 1 < n < N, 
are bounded for sufficiently large N. This can be easily 
inspected by numerical integration of the Hill's equation, 
taking into account that the boundedness of X^ (9) is 
equivalent to the condition |X+ (T/2)| < 1. We restrict 
ourselves to the summation from to N in (fT6|) . In de- 
ciding on N, it is necessary to take into account that the 
related error in Q must not exceed 0(>t). Below we take 
k = 0.005, N = 9, and set X (0) = X Of) {0) = for 
simplicity. 

Figure [3] shows the examples of the admissible defor- 
mations calculated for three different values of a max . 
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FIG. 3. Admissible deformations of the initial pulson's profile 
calculated by formula (16]) with X {0) = and X„(0) p3)l. 

In Figs. 0]and[5]we compare our solution ([9j)-(fT3"|) (solid 
lines) with the results of direct numerical integration of 
the EKG system (indicated by dots). We started with 
one of the admissible deformations of the pulson's profile 
that we have found (see Fig. [3]). Oscillations of the scalar 
field and metric at the center of the pulson are shown in 
Fig. |4j Figure [5] shows the pulson's and metric's profiles 
taken in some intermediate moment of time. 

We have performed the Fourier analysis of the scalar 
field oscillations obtained by numerical integration of the 
EKG system. The resulting spectrum shown in Fig. EJa) 
demonstrates periodicity with high accuracy. 

Then we violated the condition (f3"3"f by tripling X\ (0) 
that was calculated before, and integrated the EKG sys- 
tem again. As expected, the resulting field oscillations 
were found to be non-periodic. The corresponding spec- 
trum is presented in Fig. H]Jb). Nonperiodicity manifests 



itself as additional peaks in the spectrum which are ab- 
sent in Fig. [nja). 

To clarify the meaning of gravity, we used the obtained 
initial conditions (I36p . (|3T[) with x as a formal parameter 
for the numerical integration of the nonlinear KG equa- 
tion ([5]). The solution was found to be non-periodic, as is 
clear from its spectrum which is shown in Fig. [7] We thus 
conclude that it is because of gravity that the periodic 
pulsons of the considered non-Gaussian shapes exist. 



IV. CONCLUDING REMARKS 

Thus we have demonstrated the existence of long-lived 
time-periodic pulsons in the EKG system. These pulsons 
differ from the non-gravitating ones in their shapes and 
frequencies and exist only due to gravitational effects. 

The question arises as to whether these gravipulsons 
are stable. While the stability analysis is out of the scope 
of the present work, it is worth noting that the stability 
of the solution (|9]), and hence (TTTT) . is determined by the 
stability of the solutions of the non-homogeneous Hill's 
equation X n (9) involved in (|16[) . In turn, as it follows 
from (|M|) . the stability of the general solution X n {9) is 
determined by the behavior of the functions X£(9) and 
X-{8). 

It is clear that all solutions X n (9) satisfying the initial 
conditions ([33]). (f34|) are unstable in the resonance case 
fi n > 0. Indeed, any perturbation of the initial values 
X n (0), X ng (0), determined by J33]), flM|), leads to the 
appearance of terms ~ exp(^„#) on the right-hand side 
of (|24|) . thus making the corresponding function X n (9), 
and hence the solution @, exponentially growing in time. 

On the other hand, in the non- resonance case fi n = iv ni 
the functions X+(9), X~{9) as well as X n (8) in (JH) are 
bounded, and a small perturbation of the initial condi- 
tions (|3"3")l . results in the appearance of only small 
oscillating terms in X n {9). So we can expect that the so- 
lution © is stable if all modes X n {9) are non-resonant. 
The question is, does any value of a max exist such that 
all modes E n = An (n > 0) are stable? 

A collection of the peaks \i n > with n = 1, 2, 10, 
shown in Fig. [51 demonstrates the existence of numer- 
ous stability intervals separating the instability ones on 
the aj^ax axis. All modes n = 1,2,..., 10 are stable in 
the gaps between the peaks, and thus /i n = iv n . How- 
ever, if we take into consideration additional modes with 
10 < n < N, supplementary peaks must be added to this 
plot. Some of the new peaks will be overlapped by the 
existing ones, but the rest will fall within the stability 
intervals and erode them. Nevertheless, narrow stability 
gaps remain visible on the abscise axis even in the case 
of large N. 

While we have no proof that some gaps of stability sur- 
vive as N goes to infinity, one should take into account 
that the amplitude of the peaks in Fig. [^decreases with 
increasing n, and in any case, narrow intervals on the 
^max axis can be found where only high-n modes are un- 
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FIG. 4. Oscillations of the scalar field (top panel) and metric coefficient (bottom panel) at the center of the pulson for x = 0.005, 

dmax = 0.64. 
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FIG. 5. Profiles of the scalar field and metric coefficients. 



FIG. 6. Fourier spectrum of <j>(t, 0) for admissible (a) and 
inadmissible (b) initial conditions. 



stable. We refer to them as intervals of quasi-stability. 
Indeed, while the solution (|9]) with a max falling in one 
of these intervals is unstable, this instability evolves very 
slowly, and the gravipulson still remains a long-lived ob- 
ject. Moreover, as it was demonstrated in our simulation 
|20| , in the case of the non-gravitating pulson, the nonlin- 
ear stage of instability saturates very quickly, resulting in 
a slightly modified pulson which remains a compact os- 
cillating object. We expect the same instability behavior 



in the case of gravipulsons also, at least at small tt, even 
if this instability is caused by the action of some other 
perturbative objects around them. 

A few words about possible astrophysical applications 
of the obtained solution are in order. There are a number 
of papers where scalar solitons are considered as models 
of galactic halos in hopes of explaining the observational 
flatness of the rotation curves (see., e.g., [3l[ and refer- 
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FIG. 7. Fourier spectrum of <j>(t, 0) obtained from the solution 
of the nonlinear KG equation (JSJ with the initial conditions 
(|36|l . (|37[) . All parameters are the same as in Figs. (J4j) , ([5j . 



ences therein) . It is easy to see that then the energy den- 
sity of a scalar field must not decay faster than r~ 2 . Ev- 
idently, our solution does not satisfy this criterion. How- 
ever, if a galactic halo is not a single soliton-like object, 
but is an ensemble of dark matter lumps, of so-called 
" mini-MACHOs" [32j], the gravipulsons may be reason- 
able candidates for these compact constituents. In this 
case the gravipulson masses (1551) need to be limited by 
the condition M < 1O _7 M0 following from microlensing 
data [33[ . This constrains the amplitude of the gravipul- 
sons and the parameters of the potential © . 
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